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Abstract 

A systematic method for constructing increasingly precise pulse sequences to refocus spin evo- 
lution in the presence of large resonance offset is presented. Error cancellation relies, in part, 
on the spinor transformation property that yields a phase (—1)^ for a rotation of 2fc7r (integer 
k). The sequences only require the ability to perform vr pulses about two opposite axes, lending 
simplicity to experimental implementation. Very high fidelities with the identity operator can be 
achieved for resonance field offsets comparable to the resonant field amplitude, particularly when 
using appropriate composite vr rotations as base pulses. The sequences possess good robustness 
to systematic amplitude errors, and for certain choice of base pulse, a systematic under-rotation 
allows the bandwidth of high-fidelity refocussing to be extended to offsets larger than resonant 
field amplitude. 

PACS numbers: 03.67.Pp, 76.30.-v, 76.60.Lz 
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I. INTRODUCTION 



Utilizing spins as quantum bits (qubits) for solid-state quantum information processing 
(QIP) Q poses many practical challenges, but particularly acute are those concerning the 
manipulation and coherence of spin states given limited controls and a real environment. 
The problem of controlling a single spin with an unknown, static resonance field offset is 
often encountered in magnetic resonance, and is pertinent to recent experiments on single 
electron spin resonance (SESR) in semiconductor quantum dots containing lattice nuclear 
spins {3,4,5]. In this paper, we shall focus on the special problem of refocussing spin evolu- 
tion (i.e. generating a pulse sequence to approximate the identity operator) under resonance 
field offsets comparable to the control field amplitude. To date, SESR experiments using 

n nn 

applied microwave fields [3'] or electric-dipole spin resonance [4, |5[ have operated in a regime 
in which the distribution of random, quasi-static nuclear fields gives rise to unknown 
resonance frequency offsets larger than the available Rabi frequencies. Recent experiments, 
however, suggest that technological improvements could make significantly faster spin rota- 
tions available j^, so that control amplitude and offset may be comparable. The ability to 
refocus spin evolution with high fidelities in this setting would provide a means for precise 
measurement of decoherence times and for preserving qubit coherence over measurement en- 
sembles. A robust method that is simple to implement experimentally would be especially 
desirable. 

One approach to this problem is to focus on improving the fidelity of a single vr rotation 
by designing composite pulses that cancel errors to successively higher orders, at the expense 
of exponential growth in pulse length with respect to leading-order error 6|, |7| . The optimal 
pulse sequence for a desired period T would then heP — r — P — r, where 2(r + tp) = T and 
P represents a composite vr-pulse of duration tp. Optimal performance would be expected 
for tp = T/2, i.e. use of the highest fidelity (longest) composite pulse. Conceptually this 
is an attractive approach, but so far it has proven difficult to construct very high order 
composite pulses to refocus off-resonance errors , as the finite excitation bandwidth 
for a given resonant field amplitude (and control bandwidth) ultimately limits the achiev- 
able accuracy for arbitrary offset. While these fundamental limits apply to any approach, 
other methods have been developed to generate robust decoupling sequences from limited 
controls, such as Eulerian cycles on a group and concatenated dynamical decoupling 
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sequences with imperfect pulses 10|. This paper demonstrates an alternative method of 
approximating the identity operator in the special case of large resonance offsets, that al- 
lows errors to be cancelled to arbitrary order in principle, and that should lend itself well to 
experimental implementation. We assume a base 7r-pulse P with some fixed fidelity profile 
versus resonance offset, and then construct sequences using P and its phase- reversed twin 
pulse P (i.e. a rotation about the opposite axis in the Bloch sphere). Similar to the case of 
a single composite pulse, it is evident that a constant amplitude (windowless) sequence of 
pulses can make maximum use of the available excitation and control bandwidths. Given 
a 27T pulse duration 1/z/i (z/i is the resonant field amplitude) that is an integral fraction 
of T, we look for sequences of = 2z/iT vr-pulses that will best approximate the identity 
operator. Our algorithm for constructing these error-compensating sequences relies explic- 
itly on the transformation property of spinors that an ideal rotation of 2k7i generates an 
evolution operator (— 1)'^1, where 1 is the identity matrix and > is an integer. As longer 
sequences are constructed, higher order errors due to resonance offset are canceled, so that 
extremely high fidelities are reached for e << 1, and refocussing bandwidths can be pushed 
out to e ~ 1, limited ultimately by the bandwidth of control (a property of the chosen base 
pulse P). Since the sequences are built up from only two base pulses, P and P, experimental 
implementation is simple compared to high-order composite pulses or numerically optimized 
control sequences. We also show that these sequences can be applied in windowed mode 
(i.e. with interspersed, equal delays) without appreciably degrading the fidelity. A benefit 
of multi-pulse refocussing relevant to QIP applications is suppression of homogeneous de- 
phasing due to spectral diffusion of a coupled bath, when the bath correlation time is much 
longer than the spin nutation period. Indeed, this is the expected working regime for an 



electron spin qubit coupled to a nuclear spin bath in a quantum dot 



ll|. However, we set 



aside the detailed analysis of these sequences under dynamic errors for future work. 



A. Brief review of composite pulses 

Treating resonance offset errors in the context of average Hamiltonian theory HQ , we 
can write the time- dependent propagator of a pulse as U{t) = Uideai{'T)Uerror{'T) where r is 
the pulse length and Uerror{t) = Te'^-^o ^(*')'^*'. The interaction frame error Hamiltonian is 
V{t) = U~^l^i(t){neh'i(Tz)Uideai(t), where e is the fractional resonance offset with respect to 



3 



resonant field amplitude z/i, and az the z Pauli matrix. Tycko [13| has designed wideband 
composite pulses by numerically minimizing the lowest order average Hamiltonian terms 
in the Magnus expansion for Uerrorij), a conceptually elegant but computationally tedious 
method. Other methods have been used to construct composite rotations robust to both 
off-resonant and control field amplitude errors, such as the construction of broadband pulses 

n n r 

by Wimperis and the quaternion-based formalism of Cummins and Jones [15|, |16 |. 
i^'undamental concepts and properties of composite pulses are laid out in a review by Levitt 
Gj. Brown et al. devised formal methods for generating arbitrarily accurate sequences 
(for arbitrary single qubit rotations) that compensate for general errors j?], and provided 
specific constructions for pulse length errors but not for off-resonant errors. More recent 
work on implementation of robust state transformations and unitary operators has been 
carried out in the context of optimal control theory, e.g. use of the GRAPE method 17|, 



18l . Il9| for numerical optimization of pulse waveforms. Given some set of pulses with fixed 
precision, we are faced with the challenge of how to design sequences of many pulses that 
will fully utilize the control bandwidth so that errors do not accumulate. Indeed, an optimal 
refocussing sequence should always be non-repeating, and tailored to the desired refocussing 
period. The sequences derived below use symmetry and the spinor transformation property 
to strategically cancel errors as longer sequences are constructed. 



II. SEQUENCES 
A. Definitions 

In the following we will represent 2x2 unitary matrices as f/ = al+ih-a^ where {a^, ay, az) 
are the Pauli matrices and + |6p = 1. Turning on a resonant field of strength uji along 
the X direction produces a unitary: 

= a{e, e)l + iK{e, e)a, + ihz{9, e)(7, (1) 

where eui is the resonance offset, hz ~ 0(e) and ~ sin(6'/2) + C(e^). We denote the 
unitary corresponding to a phase-reversed resonant field as f/j. = al — ihxCJx + ibzOz- We 
shall be interested in general unitaries of two forms, approximate rotations by vr and 2k% 
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(A; = 0, 1, 2, . . . ), respectively: 



TT 



b 



a 



b 



a 




(2) 



(3) 



where the functions denoted g, g' are the error terms that, in general, can be at any order 
in e (or e^). The rotation axis for C/^ is xcos0 + y sin0 for arbitrary 0. 

Observe that the combined unitary UaUa, which we term antisymmetric, has leading- 
order az coefficient b'^ = 2abz, and transverse coefficients by = 2bxbz, b'^ = 2bybz. The 
az term is the same as for the symmetric version UaUa, whereas the transverse terms are 
smaller by the factor bz- Hence, the transverse error for an approximate identity sequence 
is significantly reduced by using the antisymmetric structure, but only in proportion to the 
existing cr^ error. Therefore the total error of an identity sequence ultimately depends on 
how far the o"^ error can be reduced by construction of the sequence; it can be guaranteed 
that the cXz error will be dominant for e << 1 at the expense of doubling the sequence length. 
We also note that a symmetric- antisymmetric (SA) sequence of the form UaUaUaUa reduces 
the transverse error by an additional factor bz so that it is ~ 0{bx,ybl). 

B. Construction algorithm 

1. Find two sequences 5*1, 5*2 that have equal leading order error: 5*1 to generate a 2kTT 
{k even) rotation and 5*2 to generate a 2A;'7r {k' odd) rotation; 

2. Note that S^Si Ri 1 + ... + 2^z'^z and S2S2 f« 1 + ... — 2{-,(T^, where ^2 is the leading 

order error of both Si and 6*2, and the — in the second expression is due to the 
spinor rotation property; then form the combined sequence A = S1S1S2S2 which has 
a reduced az error ~ 0(e^) due to cancellation of ^z ~ ^(£^) terms, where p > q. 

3. From A we can form two useful sequences, antisymmetric AA and SA sequence AAAA. 

4. For the next larger sequence, note that ^4 is a Ott rotation, and find a sequence B of the 
same length that generates a 2k'7r {k' odd) rotation and has the same leading-order 
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Sequence 


0{5^) 


0{8y) 


C('5z) 


kmax ) / = 0.99 {simple, 3-pulse, 7-pulse} 


'2'=PP 


hxO. 


a ~ 




0.25 , 0.47 , 0.72 


%'=PPPP 


bxOpbz 




abz ~ 


0.32 , 0.54 , 0.74 


'8'=PPPPPPPP 


a%l 


a%l ~ 


abz ~ 


0.25 , 0.49 , 0.72 


'W=AA 


aHi 






0.43 , 0.77 , 0.79 


'32'=AAAA 


aHl 






0.36 , 0.60 , 0.76 


'%A'=CC 


a%l 






0.73 , 0.80 , 0.87 


'128'=CC'C'C 


a%f 




abl ~ e'' 


0.45 , 0.79 , 0.86 


'256'=FF 






abl ~ 


0.72 , 0.79 , 0.86 



TABLE I: Characteristics of successively larger 'optimized' approximate identity sequences, where 
'N' indicates the number of base %(7ry) pulses P{P)- Residual orders of error terms 6 for (Tx^ 
and az operators are given in the notation of equation 2. The maximum fractional offsets 1 6772^0; | 
for 10~^ infidelity with the identity operator are shown using a simple vr-pulse, a 3-pulse composite 



TT [6t] , and a 7-pulse composite tt 
section IIIBI 



13l | as base pulses. Sequences ^4, C and F are defined in the text, 
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FIG. 1: Numerical results for logarithmic infidelity versus fractional resonance offset e for the 
(windowless) pulse sequences listed in Table 1. e = 1 corresponds to the resonance offset equal 
to the control field amplitude. Horizontal dotted lines indicate fidelities of 0.99 and 0.9999. Since 
these are windowless sequences, the total time for each sequence is the base pulse duration tp x N, 
where is the number of base pulses, (a) Base pulse P is a simple vr pulse, (b) P is the three-pulse 

n n 

composite vr rotation {tt /2)x{3ir/2)y{TT/2)x [6]. (c) P is a 7-pulse composite vr rotation |l3l |. 

az error. (For k' = 1, B can be found that only differs from A by the phase reversal 
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of one particular base pulse); 
5. label 5'i = ^4 and S2 = B, jump to step 2. 

For example, labelling P = U„ and P = Uj^, the two elementary sequences PP and PP can 
be used for step 1: 

Si = PP pa 1 - i2bzax + i2hxhz(Ty + i2ahz(Tz (4) 
5*2 = PP — 1 + i2ahx(Tx + i2aay + i2abz(Jz (5) 

Note we have taken by = 1 + C(e^) (i.e. P is an approximate tt pulse along cr^), without loss 
of generality. The leading coefficient of 1 in the PP sequence is —1 owing to transformation 
of a spinor under 27r rotation. We have: 

SiSi 1 + ... + iAabzaz (6) 

8282^1 + ...-i4:abzaz (7) 

so that 

A = S1S1S2S2 1 + ... + 0{a% + abDuz (8) 

The 16-pulse antisymmetric sequence AA, where A = S1S1S2S2, has cr^ error ~ 0{a%z + 
a6^) ~ and transverse (cr^, cTj^) errors ~ 0{a?b^^) ~ (^(e''). 

Next, we take A as the Ott rotation and find another sequence B for a 27r rotation with 
the same cr^ error. It can be verified that the sequence PPPPPPPP (i.e. the same as A 
but with the first pulse reversed in phase) satisfies this condition. We can then construct 
C = AABB which has leading-order error reduced to ~ (9(a6^) ~ 0{e'^)^ yielding the 
useful 64- and 128-pulse sequences CC and CCCC . The leading transverse errors of CC* are 
~ 0{a%^^) ~ 0{e^^). For the next iteration, a 27r pulse that complements C can be found, 
for example, D — AAAB, yielding a 256-pulse sequence FF where F — CCDD. Hence, 
the process can be iterated to any final sequence length of size 2^+" base pulses, where n 
is a positive integer. Since error cannot be reduced for sequences smaller than 8 pulses 
(relative to the elementary sequence PP), the 'optimal' sequences for 4 and 8 pulses are 
simply the antisymmetric and SA sequences PPPP and PPPPPPPP , as listed in Table 
1. 
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FIG. 2: Numerical results for windowed sequences, using the 3-pulse composite vr rotation, (a) 
Each base pulses is preceded by a delay period At = S/vi, where vi is the control field amplitude 
in Hz. (b) At = (c) At = Delays give rise to small amplitude oscillations of fidelity 

versus offset, but do not appreciably degrade the overall sequence performance, for any delay. 

C. Performance: numerical simulation 

Figure [1] displays numerical results for the sequence fidelity with respect to the identity 
operator, defined as /(e) = |Tr(f/Jl)p/4 for a sequence that generates an offset- dependent 
unitary U^. Three types of base pulses are used: a simple vr pulse, and two types of composite 
vr rotations designed to compensate off-resonance errors. These are the 3-pulse composite 
rotation (7r/2)a;(37r/2)j^(7r/2)i,. due to Levitt [6], and the 7-pulse composite rotation due 



to Tycko 131] . Clearly, the overall performance of our sequences depends strongly on the 
leading-order error terms of the chosen base pulse. For both simple and composite vr pulses 
P we generally have a ~ 0{e'^) and ~ 0{e), and for PP, a' ~ 0{e^) and b'^ ~ ^{^^)- 
For the simple pulse, all coefficients are of order unity, whereas for composite pulses they 
can be much smaller. For the Tycko composite pulse, a ^ lO^^e^ and bz ~ 10~^e, and for 
PP, b'^ ~ abz ~ 10~^e^. Since each final sequence is either antisymmetric or SA, we are 
guaranteed that the residual error will be dominated by the error for e << 1. 

Figure [U shows that the range of offsets over which very high fidelities can be achieved 
improves significantly with sequence size. In particular, for the simple pulse and three-pulse 
cases, the 'bandwidth of fidelity' at the / > 0.9999 range improves progressively for every 
sequence of length 2"^, for m even. With the exception of m = 2, these are the antisymmetric 
sequences for which error is reduced compared to the previous sequence (see Table 1). 
Shght losses of fidelity are seen for the AS sequences with m odd, as they are constructed to 
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further suppress transverse error, but allow the dominant error to double. For the 7-pulse 
Tycko sequence, the off-resonance compensation of the base pulse is particularly good, and 
all the non-trival sequences (m > 2) result in infidelities ~ O{10^^) to at least |e| ~ 0.67. 
The improvements with sequence size at large offsets are more modest for this base pulse, for 
reasons discussed below. On the other hand, the longer sequences (m > 4) yield substantial 
improvements in precision (to error ~ O{10^^ — 10"^)) for offsets up to |e| ~ 0.6 — 0.7. 

For infidelity < 10^^, the optimal sequence length for all three base pulse types is the 
64-pulse sequence; from Table 1, we see that maximum offsets \emax\ of 0.73, 0.80, and 
0.87 can be tolerated for the simple, 3-pulse, and 7-pulse base pulses, respectively. For 
larger sequences, although error orders can continue to be reduced, the coefficients of the 



leading-order terms grow linearly with sequence size[2l| so that increasingly small gains in 
\^max\ ~ 1 are obtained. From Table 1, we see that the order of az error is 0{e"^) for 
sequences of length 2"^ (m odd), i.e. a logarithmic increase with sequence size. Transverse 
error orders also increase logarithmically, but at a much faster rate ((9(e^™~^) for m odd). 
However, the linear increase in magnitude of error coefficients severely limit the gains at 
large |e| for sequences larger than an optimal size (~ 64 pulses). It is reasonable to expect 
these limitations from the perspective of control bandwidth; extending refocusing to larger 
offsets would necessarily require either larger pulse amplitude (larger excitation bandwidth) 
or larger control bandwidth, i.e. sequences that shift pulse phases on shorter timsecales. 
For the sequences derived here, control bandwidth is independent of sequence size, so the 
effective bandwidth of fidelity saturates. Control bandwidth can be increased by choice of 
base pulse; both composite pulse types used here require larger control bandwidth than a 
simple pulse, and thus produce significantly better performance. 

Figure [2] shows fidelities versus offset for windowed sequences, i.e. sequences with an 
added free evolution delay period adjacent to each base pulse. Results are shown for the 
3-pulse composite base pulse, for three different delay periods. These free evolution periods 
cause small-amplitude fidelity oscillations as a function of offset, but do not degrade the 
overall sequence performance appreciably, regardless of the delay time chosen. This can 
be seen by comparison with the windowless sequences, figure [TJo. Therefore, any of the 
sequences can be applied for arbitrary times T > Ntp {tp is the duration of P, and 
sequence size) with results similar to the corresponding windowless case. 

Another property of interest for these sequences is the robustness to systematic pulse 
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FIG. 3: Robustness of sequences versus systematic amplitude errors, for the 3-pulse composite 
base pulse. Sequence fidelity is displayed as a contour map (grayscale), with / = 99% as the cutoff 
(white areas represent / < 99%). The horizontal axes are amplitude scale factors, and the vertical 
axes are fractional resonance offset e. It can be seen that robustness to amplitude errors generally 
increases with sequence length, in a similar fashion to offset robustness. 

amplitude error. For zero resonance offset, the sequences are perfectly robust since they 
correspond to rotations of Ovr. Hence, we expect reasonably good amplitude robustness in 
the range that the pulses continue to compensate well for offset errors. This expectation is 
borne out in numerical simulations: figure [3] shows sequence fidelity versus both resonance 
offset and an amplitude scale factor, for several sequences using the 3-pulse composite base 
pulse. Robustness to amplitude errors increases with sequence size in a similar fashion to 
offset robustness. For the 3-pulse composite vr, the patterns demonstrate that low-amplitude 
errors, or under-rotations, are more forgiving than over-rotations. In this case, given a sym- 
metric distribution of amplitudes, the optimal average amplitude for maximizing fidelities 
over the distribution is ~ O.dui. However, this behavior is dependent on choice of base pulse; 
for a simple pulse, the fidelity versus amplitude is more symmetric and centered on scale 
factor 1. This underrotation property is also shared by the Tycko (7-pulse) composite pulse, 
and can be used to increase the bandwidth of fidelity at the expense of reducing fidelities 
for small offsets. Figure Hti shows the fidelity of the 64-pulse sequence using the Tycko pulse 
versus offset and amplitude scaling. The vertical dotted line at scale factor = 1 indicates 
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the cut shown in the figure [T}:;, i.e. 'correct' amphtude. The left dotted fine at scale factor 
^ 0.8 indicates a wider region of offset over which sequence fidelities / > 0.999 are obtained. 
Hence, by underrotating each base pulse P by ^ 0.8, a wider bandwidth of refocusing can 
be achieved in a certain fidelity range. This is somewhat surprising, since lower amplitude 
corresponds to a smaller bandwidth of excitation. Figure shows that the same holds true 
for underrotating by scaling down the timesteps of the base pulse while maintaining the 
original pulse amplitude. Infact, this method allows even wider bandwidth refocusing: for 
the 64-pulse sequence, fidelities / > 0.999 are achieved across the range < |e| < 1.2 when 
scaling the timesteps by 0.8. This is can be understood as an increase in the bandwidth of 
control. The tradeoff is larger infidelity at small offsets, i.e. / ~ 0(10"'^) for the 64-pulse 
sequence. 




1 1.5 ■ 0.5 1 

amplitude scale factor fractional resonance offset e 



FIG. 4: Larger bandwidth refocusing can be achieved for the Tycko (7-pulse) composite base pulse 
by underrotating, at the expense of reduced fidelities in the low-e regime, (a) For the 64-pulse 
sequence, by reducing the pulse amplitude to ^ O.Si^i (left dotted line), the threshold |emax| for 
fidelities / > 0.99 can be pushed closer to 1. (b) Similarly, by maintaining the original pulse 
amplitude but rescaling the pulse timesteps by 0.8, fidelities / > 0.999 can be achieved to |e| 1.2, 
with the tradeoff of flatter performance for < |e| < 1.2 compared to the original timestep scaling. 



III. CONCLUSIONS 

In conclusion, we have shown a new method for constructing error- compensating se- 
quences to approximate the identity operator for the case of resonant field strength com- 
parable to resonance offset. The sequences require only one base pulse, together with its 
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phase reversed twin, lending simplicity for experimental implementation. The construction 
algorithm makes explicit use of the spinor transformation property as a tool for canceling 
offset errors. To construct the sequences, we first reduce errors using the spinor rotation 
property, then reduce transverse error (in proportion to o"^ error) by anti-symmetrization, 
repeating iteratively. The process can, in principle, be carried out to any error order, but 
the logarithmic elimination of error orders relative to growth of error coefficients limits the 
gains in refocussing bandwidth for very long sequences. 

We expect these sequences to be useful for maintaining coherence in any qubit devices for 
which there exist wide, static distributions of qubit splitting and limited control field ampli- 
tude, or for which many measurements must be averaged in the presence of slowly drifting 
resonance frequencies. One particular example is the presence of random (quasi-static) nu- 
clear fields acting on electron spin qubits in quantum dots; such wideband refocusing with 
limited controls could provide a useful means to effectively decouple the nuclear system from 
the qubit. It will be interesting to investigate the efficiency with which slow evolution of 
the nuclear fields, driven by spectral diffusion, spoils the error cancellation properties of 
these sequences, i.e. whether such sequences may provide any benefits for decoupling from 
a dynamic environment. 
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